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摘要 : 新 能 源 多 馈 入 系统 存在 由 于 锁 相 环 带 宽 与 电网 强度 不 匹配 导致 的 次 / 超 同 步 振荡 问题 ， 
度 以 保障 系统 稳定 运行 。 为 此 ， 
出 力 或 端口 电压 偏离 额定 值 的 影响 , 提出 适用 于 非 额定 工 况 下 新 能 
并 提出 了 广义 运行 短路 比 指标 。 本 文 首 先 分 析 了 新 能 
量 的 相互 作用 关系 ; HK, 


要 准确 量化 电网 强度 


规律 ， 近 似 解 耦 设备 动态 特性 与 端口 电气 
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BUR AAR IRIE AR PIS 
EE 源 系 统 电网 强度 的 量化 方法 ， 
E 源 系统 运行 方式 变化 时 设备 的 外 特性 变化 
研究 了 考虑 电气 量变 化 时 


该 文 考虑 新 能 


的 网 络 侧 导 纳 矩 阵 的 可 解 耦 性 ， 提 出 了 考虑 运行 方式 的 广义 运行 短路 比 指标 和 相应 实用 化 计算 
方法 ， 据 此 拓展 了 广义 短路 比 理论 的 应 用 范围 。 最 后 ， 基 于 算 例 说 明 所 提 方 法 的 有 效 性 。 
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引言 

高 比例 新 能 源 、 高 比例 电力 电子 设备 接 入 电网 
导致 “ 双 高 系统 ”强度 相对 变 弱 [61]， 可 能 引起 由 
于 锁 相 环 带宽 与 电网 强度 不 匹配 导致 的 次 / 超 同 步 
振荡 问题 时 品 。 此 外 ， 由 于 新 能 源 实际 出 力 的 随机 
性 加 和 交流 电网 开 停机 或 故障 等 导致 的 系统 网 络 

架构 改变 带 来 的 电压 扰动 , 使 得 实际 系统 中 新 能 源 
并 网 端口 特性 复杂 多 变 ,影响 系统 的 小 干扰 同步 稳 
性 。 

为 保证 电力 系统 安全 稳定 运行 , 准确 度量 多 馈 
入 系统 的 电网 强度 和 稳定 裕 度 至 关 重 要 。 目 前 工程 
上 主要 采用 基于 短路 比 (Short Circuit Ratio, SCR) 
的 方法 TL。 根据 不 同 适用 场景 ， 多 馈 入 短路 比 
可 以 分 为 两 大 类 。 
第 一 类 指标 适用 于 系统 规划 阶段 ， 如 CIGRE 
组 织 基于 母线 等 值 思想 ,通过 引入 多 直流 交互 因子 
(multi-infeed interaction factor, MIIF) 提 出 的 多 馈 入 
短路 比 (multi-infeed short circuit ratio, MISCR) 指 标 
D2709， 该 指标 定义 简单 ， 缺 乏 严 格 的 理论 基础 ; 
文献 [15]-[16] 基 于 模 态 解 耦 思想 提出 的 广义 短路 
tk (generalized Short Circuit Ratio, gSCR ) 和 临界 
广义 短路 比 (Critical generalized Short Circuit Ratio) 
指标 通过 严格 的 数学 变换 , 将 符合 条 件 的 多 馈 入 系 
统 近 似 解 看 为 若干 个 单 馈 入 系统 ,从 而 将 多 馈 入 系 
统 的 稳定 性 分 析 聚 焦 于 最 弱 单 馈 入 系统 , 便于 将 单 
馈 入 系统 中 短路 比 的 经 验 移 至 于 多 人 馈 入 系统 , 应 用 

上 分 方便 。 
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次 / 超 同 步 振 荡 ; 非 额 定 工 况 ; 广义 运行 短路 比 ; 可 解 夺 性 


另 一 类 适用 于 实际 运行 阶段 , 它 考虑 了 系统 非 
额定 运行 方式 对 电网 强度 的 影响 ， 如 文献 [17] 提 出 
的 适用 于 MMC 系统 的 临界 运行 短路 比 (Critical 
Operating Short Circuit Ratio，COSCR) 指 标 ， 将 交 
流 母 线 电 压 约 束 条 件 下 最 小 的 临界 短路 比 作为 
COSCR， 人 然而 在 实际 系统 存在 电压 波动 时 ， 利 用 
该 指标 进行 稳定 性 分 析 可 能 出 现 误 判 的 情况 ; 文献 
[18] 提 出 了 适用 于 多 直流 馈 入 系统 静态 电压 稳定 
性 分 析 的 等 效 运行 短 路 比 (equivalent operating 
short circuit ratio，EOSCR) 指 标 ， 但 该 指标 未 考虑 
电力 电子 设备 间 的 相互 作用 ， 评 估 可 能 过 于 乐观 ; 
文献 [19] 提 出 的 等 价 多 馈 入 有 效 短路 比 指标 
(equivalent multi-infeed effective short circuit ratio, 
EMESCR) 综 合 考虑 了 阻抗 、 多 馈 入 相互 作用 因子 、 
无 功 功 率 和 节点 电压 变化 等 实际 因素 对 系统 稳定 
性 的 影响 , 但 该 指标 定义 复杂 , 且 未 提出 对 应 临界 
指标 ， 难 以 量化 系统 稳定 裕 度 。 

为 准确 度量 系统 在 实际 运行 阶段 的 电网 强度 
和 稳定 裕 度 ， 和 针对 锁 相 环 主导 的 小 干扰 稳定 问题 ， 
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本 文 拓展 了 经 过 严格 数学 推导 的 广义 短路 比 理论 
的 应 用 范围 ， 提 出 广义 运行 短路 比 (generalized 


Operation Short Circuit Ratio, gOSCR) 和 临界 广义 
运行 短路 比 (Critical generalized Operation Short 
Circuit Ratio, CEOSCR), 近似 解 耦 设备 动态 特性 与 
端口 电气 量 的 相互 作用 关系 , 并 基于 摄 动 理论 研究 
了 非 额 定 运 行 方式 下 电网 强度 的 量化 机 理 、 方 法 ， 
提出 了 指标 的 解析 和 实用 化 计算 方法 ,实现 了 多 馈 
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入 系统 实际 运行 阶段 电网 强度 和 稳定 裕 度 的 准确 
量化 评估 。 最 后 ， 分 别 基 于 MATLAB/Simulink 平 
台 的 仿真 模型 和 CloudPSS 平台 的 哈密 新 能 源 场 站 
等 比例 缩小 电磁 仿真 模型 说 明了 所 提 指 标的 有 效 
性 。 
1 单 馈 入 系统 运行 阶段 小 干扰 稳定 性 
1.1 单 馈 入 系统 短路 比 

考虑 变 流 器 接 入 一 个 无 穷 大 系统 , 其 控制 框图 
见 附录 A 图 A1。 目前 新 能 源 小 干扰 稳定 问题 主要 
为 变 流 器 锁 相 环 (phase locked loop，PLL) 控 制 环节 
主导 的 次 / 超 同步 振荡 ， 其 频段 大 致 在 中 低频 段 
(10~100Hz) 20， 此 外 ， 根 据 《 光 伏 发 电站 接 入 电力 
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= U(H, +sL,)(1+UH,,) 


下 面 分 析 传 输 网 络 的 动态 。 
交流 电网 的 动态 模型 为 (端口 电流 以 流出 节点 


系统 技术 规定 》(GB/T 19964-2012)、《 风 电场 接 入 
电力 系统 技术 规定 》(GB/T19963) 中 分 别 对 其 无 功 
容量 的 要 求 ， 系 统 稳定 运 行 时 ， 风 电 、 光 伏 等 新 能 
源 功 率 因 数 约 为 1。 因 此， 基于 文献 [20]-[21] 中 对 
变 流 器 内 部 控制 系统 多 时 间 尺 度 特 性 的 描述 ,在 中 
低频 段 下 , 变 流 器 采用 开关 平均 模型 ， 主 要 考虑 电 
流 内 环 、 锁 相 环 和 传输 网 络 的 动态 , 功率 因数 为 1。 
下 面 分 析 新 能 源 设备 端口 动态 外 特性 。 
在 全 局 同步 坐标 系 下 , 得 到 变 流 器 在 中 低频 段 
下 以 端口 电压 增 量 和 端口 注入 电流 增 量 构成 的 导 
纳 传递 函数 矩阵 动态 模型 : (详细 推导 过 程 见 文献 


[20]): 
Al, — You Yon AU, (1) 
Al, Yoo Yn | AU,, 
SRS ey 
Wse (s) 
式 中 : 
_-(1-G) 
o (H, +sL,) 
Yo = Yoo =0 


U0 Gye) + PH „H, 
= U(H,+sL)(1+UH y) 

其 中 , 到 五 分 别 为 设备 输出 电流 在 全 局 同步 坐标 系 
下 的 分 量 ; Un U IDANRA (point of common 
coupling, PCC) 电压 在 全 局 同步 坐标 系 下 的 分 量 ; 
Yrsc(s) 为 变 流 器 导 纳 传递 函数 矩阵 ; Grr=1/(1+TFFs)» 
为 电压 前 馈 低 通 滤波 传递 函数 ; s 表示 拉 普 拉 斯 算 
F; HKytKa/s，Kip、Ki 分 别 为 内 环 传递 函数 的 
PI 控制 参数 ， Hp(Kpip+tKpui/s)y s> Kpips Kpui 3y 
A PLL 的 PI 控制 参数 ; U 为 并 网 点 电压 实际 值 ; P 
为 设备 实际 输出 有 功 功 率 。 

在 中 低频 段 下 ,可 忽略 电压 前 馈 动 态 PC0， 式 可 
转换 为 : 
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为 正方 向 ): 
Al 1 AU 
区 ear | 6) 
Yrer (8) 
式 中 : 
Z=Q@L, 
ra- 
-a(s) p(s) 


HA, Ynes) ABOU FE, PA Sig SRSA FE BES Z 为 交流 电 
网 的 等 效 阻抗 ， 由 为 同步 旋转 速度 ; 
B(s) =se,/(s? +03) ， a(s) =a; (s +0) 5 

基于 多 变量 频率 域 反馈 控制 理论 C]， 建 立新 
能 源 并 网 系统 的 闭环 特征 方程 : 


0 0 1 
ll a + zro) =0 (4) 


根据 式 与 短路 比 之 间 存 在 的 解析 关系 , 在 系统 
规划 阶段 ,可 利用 传统 短路 比 指标 量化 系统 的 电网 
强度 ， 短 路 比 定义 式 为 0 上: 
S.U’? 1 1 
s ZS, Sz (6) 
式 中 : Si 为 交流 短路 容量 ，S 为 变 流 器 自身 的 额 
定 容量 。 在 规划 阶段 可 认为 忌 的 额定 值 是 1， 故 
SCR =1/S,Z o 

然而 ， 当 设备 运行 在 非 额定 运行 方式 时 ， 式 表 
明 设 备 端口 电气 量 ， 即 输出 有 功 功率 和 端口 电压 会 
影响 设备 的 动态 特性 ， 从 而 影响 系统 的 小 干扰 稳定 
性 。 因 此 ， 在 实际 运行 阶段 ， 仅 利用 式 定 义 的 短路 
比 指标 来 量化 系统 电网 强度 是 不 准确 的 。 
1.2 单 馈 入 系统 运行 短路 比 

为 解 看 设备 动态 特性 与 端口 电压 和 输出 有 功 
功率 的 相互 作用 关系 ， 用 运行 短路 比 (Operation 
Short Circuit Ratio, OSCR) 完 全 表征 系统 非 额定 运 
行 方式 对 系统 小 干扰 稳定 性 的 影响 ， 本 节 基 于 变 流 
器 导 纳 传递 函数 矩阵 中 Yo 的 结构 特性 ， 对 其 进行 
简化 、 变 换 。 丈 2? 分子、 分 母 同 时 乘 以 U, 提出 P/U， 


SCR = 


式 可 变换 为 : 


0 0 
Ysc(s) = = 0 H, UH yn (6) 
(H, +sL,)(1+UH,y) 


jysc (8) 
EP, P(g) 为 提取 P/U? 后 的 变 流 器 导 纳 传递 函数 
矩阵 。 
此 时 ,， 变 流 器 单机 无 穷 大 系统 的 闭环 特征 方程 
可 重新 写 为 : 


dt| Fats) + Eser- 0)) -0 (7) 


式 反映 了 系统 的 小 干扰 稳定 性 与 SCRxUYP 的 
显 式 关 系 ， 因 此 可 以 将 运行 短路 比 定义 为 : 


U? 
OSCR = — B 8 
= (8) 


其 中 ，B 为 交流 电网 的 等 效 导 纳 。 

上 式 表 明 ， 相 比 于 传统 单 馈 入 SCR, OSCR 考 
虑 了 实际 运行 阶段 中 变 流 器 的 非 额定 运行 方式 ( 实 
际 出 力 偏离 额定 功率 以 及 外 部 电网 电压 偏离 额定 
值 ) 对 系统 小 干扰 稳定 性 的 影响 ， 因 此 能 较为 准确 
地 度量 运行 阶段 下 单 馈 入 系统 的 电网 强度 。 
13 ” 单 馈 入 系统 临界 运行 短路 比 

单 馈 入 系统 稳定 运行 时 需 满足 其 OSCR 数值 
大 于 临界 运行 短路 比 (Critical Operation Short 
Circuit Ratio, COSCR)。 与 文献 [16] 对 CgSCR 的 定 
MAW, 系统 COSCR 为 单 馈 入 系统 主导 特征 值 实 
部 为 0 时 的 OSCR， 如 下 式 所 示 : 


COSCR = arg {det (Y,.c(7@) + OSCR -F (ja) )=0} (9) 


OSCR 


下 


式 中 : arg{:} 表 示 求 方程 的 根 ，w 为 控制 器 旋转 的 角 
频率 。 

图 1(a)、(b) 分 别 描绘 了 临界 稳定 下 单 馈 入 系 
统 在 不 同 端 口 电压 和 输出 有 功 功 率 的 SCR、 
COSCR 变化 曲线 。 可 以 看 出 ， 端 口 电压 和 输出 有 
功 功 率 变 化 时 ， 系 统 在 临界 稳定 下 的 SCR 相对 变 
化 较 大 , CSCR 无 法 准确 地 刻画 实际 运行 系统 的 稳 
定 边界 ， 相 较 于 SCR， 考 虑 非 额 定 运 行 方式 对 系 
统 小 干扰 稳定 性 影响 的 COSCR FB tx HF Eels) 仍 
含有 电压 项 ， 因 此 仅 随 端口 电压 的 改变 而 改变 ， 不 
受 系 统 其 它 电 气量 变化 的 影响 ， 验 证 了 COSCR 指 
标 在 刻画 系统 稳定 边界 的 有 效 性 。 
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图 1 临界 稳定 下 单 馈 入 系统 SCR、OSCR 变化 曲线 
Fig.1 SCR and OSCR change curve of SIPES under 


critical stability with various voltage or active power. 
2 多 馈 入 系统 运行 阶段 小 干扰 稳定 性 


文献 [15] 基 于 模 态 解 耦 的 思想 ， 提 出 了 适用 于 
规划 阶段 多 馈 入 系统 电网 强度 评估 的 广义 短路 比 
指标 。 类 比 该 思想 , 本 节 将 上 述 指标 推广 至 考虑 设 
备 出 力 或 端口 电压 偏离 额定 值 的 运行 场景 , 并 提出 
一 套 实 用 化 的 指标 计算 方法 。 

考虑 图 2 所 示 由 nn 台电 力 电子 设备 馈 入 的 交 
流 系统 ， 其 中 Un 0: (i=1,...,n) 为 第 i 个 电力 电子 
设备 馈 入 节点 的 电压 模 和 相 角 ，E、o 为 外 网 等 值 
电压 幅 值 和 相 角 。 
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图 2 多 馈 入 电力 系统 


Fig2 Multiple-infeed power electronic system 


图 2 所 示 多 馈 入 系统 的 闭环 特征 方程 为 : 


Yisci(s) 0 0 
det 0 es 0 
0 0 Ysc,(s) 


其 中 , BeR™ 表示 仪 包含 设备 节点 的 网 络 降 阶 导 纳 
和 矩阵; Yysci(s) 为 第 i 台电 力 电 子 设 备 的 导 纳 传递 函 
DUE ME; © 表示 Kronecker 积 。 

类 比 式 ， 分 别提 取 Yrsci(s) 的 P/U?， 以 解 粳 设 
备 动态 特性 与 端口 电气 量 ， 式 可 变换 为 : 


+B onal =0 (10) 


Yisc: 0 y? 
det}/ 0 >. 0 rdiag()B®F(S) =0 (11) 
0 0 Ysc i 


AH: Pi:、Ui 为 第 i 台电 力 电 子 设备 的 输出 有 功 功 
率 和 端口 电压 (n 台电 力 电子 设备 使 用 统一 基准 值 
进行 标 么 ); 疡 -为 第 字 台 电力 电子 设备 的 交 (9 o 
2.1 多 馈 入 系统 广义 运行 短路 比 

当 多 馈 入 系统 中 台电 力 电子 设备 动态 外 特 
PEARL CR SCTE “TA AB”), FLY, 差异 性 仅 
来 自 于 不 同 设备 端口 电压 的 差异 性 。 

按照 《光伏 发 电站 接 入 电力 系统 技术 规定 》 
(GB/T 19964-2012)、《 风 电场 接 入 电力 系统 技术 
规定 》(GB/T19963) 对 稳 态 下 系统 电压 控制 的 要 求 ， 
考虑 节点 电压 仅 在 0.9~1.1p.u 的 安全 范围 内 变化 ， 
并 且 在 变 流 器 锁 相 环 控制 主导 的 小 干扰 稳定 问题 


中 ,系统 主导 模 态 在 锁 相 环 带宽 内 ， 则 在 系统 主导 
KAS PA |UH,,|> 1, RS AL: 
IH 
ae (12) 
(1+UH,,) 


ARH, Psa 间 差 异性 较 小 ， 因 此 可 类 比 文 献 
[16] 在 系统 规划 阶段 的 加 权 思 想 构造 等 效 同 构 系 统 ， 


以 逼近 原 系统 的 稳定 性 。 等 效 同 构 系 统 的 闭环 特征 
ai, ® A(s) + daethon] =0 (13) 

A(s) = G(s)F\(s) (14) 

G(s) =>) piYsc(s) (15) 

Pii = UV (d 6) 


其 中 , 为 nn 阶 单位 阵 ; ul 和 vw 分 别 为 qiag(U?1P)B 
关于 最 小 特征 值 4 的 左右 特征 向 量 ，xw 和 vi 分 别 
为 wr 和 wv 的 第 一 个 元 素 ， 归 一 化 后 满足 p=1。 

与 广义 短路 比 的 推导 过 程 类 似 ， 对 
diag(U; / P)B 进行 特征 值 分 解 ， 式 可 看 作 n 个 单 馈 
入 系统 特征 方程 的 积 ， 即 51; 


TIasfAaw) +41}=0 (17) 


式 中 : AA diag(U?/ P)B 的 特征 值 ，0 <4<…<4。 

引 理 1: 在 多 馈 入 系统 中 , U?/ p>0, diag(U?/ P) 
和 diag(U?1P)” 是 正定 矩阵 ， 同 时 ，B 是 对 角 元 均 
为 正 数 的 Hermitian #EpE, Æ EEE. AE, 
据 文 献 [22][26] 对 sze 特征 值 的 分 析 可 知 ， 
diag(U? / P)B 是 可 对 角 化 的 ， 且 特征 值 都 是 正 的 。 

式 表明 ， 多 馈 入 系统 的 稳定 性 可 通过 分 析 解 类 
后 的 单 馈 入 系统 得 到 ， 即 多 馈 入 系统 的 电网 强度 取 
决 于 最 弱 的 等 效 单 馈 入 系统 [3。 因 此 ， 可 将 最 弱 单 
馈 入 系统 的 OSCR, E diag(U?/ P)B 的 最 小 特征 值 4 
定义 为 多 馈 入 系统 的 广义 运行 短路 比 ， 即 : 


2 
gOSCR = om 让 | (18) 


式 中 : Ate) 表示 求解 矩阵 特征 值 。 

本 文 侧重 考虑 运行 方式 对 系统 小 干扰 稳定 性 
的 影响 ， 故 简化 设备 建 模 ,， 采用“ 同 构 系统 ”分 析 
非 额定 运行 方式 下 电网 强度 的 量化 机 理 。 当 多 馈 入 
系统 中 台电 力 电子 设 备 动态 外 特性 存在 差异 (后 
文 简称 “ 异 构 系 统 ”)， 其 疡 ,差异 性 来 自 于 端口 电 
压 和 设备 控制 结构 、 参 数 等 ， 分 析 方 法 与 “ 同 构 系 
统 ” 类 似 ， 这 里 不 再 资 述 。 

22 多 馈 入 系统 临界 广义 运行 短路 比 

多 馈 入 系统 稳定 运行 时 ， 其 gOSCR 应 大 于 
CgOSCR 。 文 献 [16] 中 CgSCR 能 够 通过 分 析 多 馈 
入 系统 对 应 等 效 同 构 系统 来 解析 或 仿真 计算 得 到 。 
同 理 ， 当 考虑 非 额定 运行 方式 下 设备 出 力 或 电压 偏 
离 额定 值 时 ， 同 样 可 通过 分 析 等 效 同 构 系统 来 解析 
或 仿真 计算 多 馈 入 系统 的 临界 广义 运行 短路 比 ， 其 
定义 式 为 : 


CgOSCR = arg {det (AC j@)+gOSCR )= o} (19) 
gOSCR 


“ 同 构 系 统 ” 需 构造 等 效 同 构 系 统 来 逼近 原 系 
统 的 稳定 性 ， 其 所 需 diag(U?/ P)B 的 左右 特征 向 量 
与 电力 电子 设备 输出 有 功 功率 、 端 口 电压 和 系统 网 
架 结 构 相 关 ， 因 此 ， 系 统 的 CgOSCR 随 系 统 运行 
工 况 的 变化 而 变化 。 

2.2.1 ”临界 广义 运行 短路 比 解析 算法 

实际 工程 中 , 在 设备 外 特性 己 知情 况 下 , 能够 
通过 解析 计算 gOSCR 和 CgOSCR 实现 系统 电网 强 
度 和 稳定 裕 度 地 准确 度量 ,分 析 流 程 图 见 附录 A 
图 A2， 具 体 步 又 为 : 

1) 根据 网 络 参 数 和 实际 系统 中 各 电力 电子 设 
备 输出 有 功 功率 、 端 口 电压 ， 得 到 diag(U?/P)B ; 
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2) 根据 式 计算 得 到 系统 gOSCR; 

3) 根据 系统 中 各 电力 电子 设备 的 外 特性 和 式 
-构造 等 效 同 构 系 统 ; 
A) 根据 等 效 同 构 系统 的 等 效 单 馈 
用 式 解 析 计 算 CgOSCR: 

5) 计算 gOSCR 与 CgOSCR 的 差 值 ， 根 据 该 
差 值 度量 该 工 况 下 多 馈 入 系统 的 电网 强度 和 稳定 
裕 度 。 
2.2.2 ”临界 广义 运行 短路 比 仿真 计算 

当 多 馈 入 系统 中 设备 参数 未 知 或 者 解析 计算 
较 复 杂 , 无 法 利用 解析 计算 方法 分 析 系 统 稳定 裕 度 ， 
但 存在 数值 仿真 模型 时 , 也 可 利用 半 实 物 仿 真 或 者 
数字 仿真 , 即 在 仿真 系统 中 搭建 等 效 单机 系统 , 调 
节 线 路 电感 和 端口 电压 , 并 结合 设备 输出 有 功 功率 
以 计算 系统 在 不 同 端口 电压 下 的 临界 运行 短路 比 ， 
即 为 多 馈 入 系统 的 临界 广义 运行 短路 比 。 

本 节 将 详细 说 明 如 何 构 造 考 虑 节点 间 电 压 存 
在 差异 的 “ 同 构 系 统 ”的 等 效 同 构 系统 仿真 模型 ， 
以 通过 仿真 计算 原 “ 同 构 系 统 ” 的 临界 广义 运行 短 
路 比 。 

已 知 变 流 器 锁 相 环 控制 环节 主导 的 次 / 超 同 步 
振荡 频段 大 致 在 中 低频 段 ， 基 于 文献 [20]、[21] 中 
对 变 流 器 内 部 控制 系统 多 时 间 尺 度 特性 的 描述 ,内 
环 控制 环节 的 响应 速度 远大 于 锁 相 环 , 则 在 中 低频 
段 内 ， 系 统 满足 : 

H UH UH 


i pll a pll 
(H,+sL,)(1+UH„) (1+UH,,) 20) 


同时 ， 根 据 Neumam AAA VE By A, GOR 
阵 4〈4ec"，C 代表 复数 域 ) 的 所 有 特征 值 满足 
4<1，(C+4 的 一 阶 近似 满足 PoP7: 


(T+4) =T-4 (21) 


系统 ， 利 


其 中 ，7 为 单位 阵 ， 维 数 与 4 相同 。 
MEAM, BE Vigo 可 简化 为 : 


ses) <|0 


ar] 
L UH 1 (22) 


[0 0 
~l0 TOE N 


| UH 


p! 


HP, Foel) 为 简化 后 的 设备 六 ce o 
内 此 ， 当 系统 为 “ 同 构 系统 ”时 ， 联 立 、 和 ， 
等 效 同 构 系统 可 变换 为 : 


A(s) =F GD) puYsals) 


0 0 
=F (s) 0 1- 1 Dn (23) 


L H pn 
[0 0 
jp- 
=F (s) 0 1- 
H nU a 


其 中 ， Uea 为 等 效 同 构 系 统 的 等 效 端口 电压 。 

因此 , 在 单机 并 网 仿真 系统 中 , 调节 线路 电感 
等 使 系统 在 设备 端口 电压 为 等 效 端口 电压 的 情况 
下 达到 临界 稳定 , 得 到 的 系统 临界 运行 短路 比 ， 即 
为 原 “ 同 构 系 统 ” 的 临界 广义 运行 短路 比 。 

综 上 所 述 ，gOSCR 和 CgOSCR 指标 考虑 了 系 
统 非 额定 运行 方式 (最 大 出 力 偏离 额定 功率 以 及 外 
部 电网 电压 偏离 额定 值 ) 对 系统 电网 强度 的 影响 ， 
因此 用 两 者 的 差 值 能 够 更 准确 地 度量 多 馈 入 系统 
的 电网 强度 和 稳定 裕 度 。 
3 算 例 分 析 及 验证 

为 验证 所 提 gOSCR 指标 的 有 效 性 ， 在 
Matlab/Simulink 平台 中 搭建 单 馈 入 和 八 馈 入 系统 
时 域 仿真 模型 , 其 中 变 流 器 采用 平均 模型 , 控制 器 
的 采样 及 控制 频率 为 10 KHz。 单 馈 入 和 八 馈 入 系 
统 如 附录 A 图 A3(a)、(b) 所 示 , 变 流 器 控制 和 网 络 
拓扑 参数 见 附录 A 表 A 1。 此外， 本 节 在 Cloudpss 
平台 上 搭建 哈密 等 比例 缩小 系统 的 电磁 仿真 模型 ， 
系统 结构 如 附录 A 图 A3 (c) 所 示 ， 变 流 器 控制 参 
数 和 网 络 拓 扑 参 数 见 附录 A K A2。 
3.1 单 馈 入 系统 时 域 分 析 

本 节 首 先 基 于 单 馈 入 系统 , 通过 仿真 分 析 阐述 
节点 电压 、 设 备 输出 有 功 功 率 变 化 对 系统 稳定 性 的 
影响 ， 为 后 文 的 理论 验证 做 铺垫 。 
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图 3 端口 电压 、 有 功 功率 振荡 曲线 

Fig.3 Terminal voltage and active power oscillation curse 
with various voltage or active power 

设置 单 馈 入 系统 中 Line=0.5p.u， 即 系统 的 短 

路 比 为 2, 已 知 临界 短路 比 为 2.17, 则 从 短路 比 的 

度 可 知 系统 是 不 稳定 的 。 然 而 ， 在 t2s 时 ， 在 


= 


外 电网 施加 0.025pu 的 电压 阶 跃 , 并 在 t=5s 时 在 外 


3.2.1 等 效 同 构 系 统 主导 模 态 近似 效果 

由 上 文 可 知 , 等 效 同 构 系 统 的 稳定 性 由 最 弱 单 
馈 入 系统 决定 ， 故 比较 等 效 同 构 系 统 的 最 弱 单 馈 入 
系统 和 原 系统 主导 模 态 之 间 的 误差 ,可 以 说 明 近 似 
效果 。 

通过 等 比例 增 大 所 有 线路 电抗 , 分 别 计算 系统 
的 gSCR 和 gOSCR 以 及 原 系统 的 主导 特征 根 轨迹 
( 记 为 “591”)、gOSCR 对 应 等 效 同 构 系 统 的 最 弱 单 
馈 入 系统 的 主导 特征 根 轨迹 〔 记 为 “s，”)、gSCR 
对 应 最 弱 等 效 单 馈 入 系统 主导 特征 根 轨迹 〈 记 为 
“s3”)， 并 分 别 计算 sa s35 si 的 相对 误差 ， 误 差 
公式 分 别 为 : 


|s.-sil sıl 
|s| 


[s -sl sı 


ô= x100% (24) 


a= x 100% (25) 


s,| 

# 1 gOSCR 对 应 等 效 同 构 系 统 的 最 弱 单 馈 入 系统 的 主 
导 特 征 根 对 比 

Table 1 Comparison of dominant characteristic roots of 


the weakest equivalent SIPES considering gOSCR 


电网 再 次 施加 0.03pu 的 电压 阶 跃 ， 设 备 的 端口 电 


gOSCR St S2 ôI 
压 振荡 曲线 如 图 3(a) 所 示 ， 单 馈 入 系统 从 发 散失 2.68 -1.05+81.25i -1.45+81.73i 0.76% 
稳 到 临界 稳定 ,最 后 到 收敛 稳定 。 由 此 可 知 ,设备 3.01 -1.60+81.67i -1.93+82.10i 0.67% 
端口 电压 变化 时 ， 利 用 SCR 指标 量化 系统 稳定 性 3.42 -2.13+82.01i -2.40+82.40i 0.58% 
和 刻画 系统 稳定 边界 存在 较 大 误差 ,可 能 会 导致 系 3.93 -2.64+82.29i -2.86 士 82.64i 0.50% 


统 稳定 性 误 判 。 
进一步 设置 Line=0.45p.u， 系 统 的 短路 比 为 
2.22, 己 知 临界 短路 比 为 2.17， 从 短路 比 的 角度 可 
知 系统 是 稳定 的 。 然 而 ， 在 t=2s 时 在 外 电网 施加 
-0.05pu 的 功率 跌落 ， 并 在 t=5s 时 在 外 电网 再 次 施 
加 -0.1pu 的 功率 跌落 , 设备 的 输出 有 功 功 率 振荡 曲 
线 如 图 3(b) 所 示 ， 单 馈 入 系统 从 发 散失 稳 到 临界 
稳定 ,最 后 到 收敛 稳定 。 由 此 可 知 , 设备 有 功 功 率 
变化 时 ， 利 用 SCR 指标 量化 系统 稳定 性 和 刻画 系 
统 稳定 边界 也 存在 较 大 误差 , 可 能 会 导致 系统 稳定 
性 误 判 。 
3.2 八 馈 入 系统 频 域 、 时 域 分 析 
系统 结构 示意 图 参考 附录 A 图 A3(b)， 


表 2 gSCR 对 应 最 弱 等 效 单 馈 入 系统 主导 特征 根 对 比 
Table 2 Comparison of dominant characteristic roots of 


the weakest equivalent SIPES considering gSCR 


gSCR U/pu 53 02 
2.92 0.9031;0.9462 -1.72+85.22i 4.96% 
3.24 0.9144;0.9527 -2.18+85.17i 4.35% 
3.65 0.9241;0.9578 -2.62+85.12i 3.83% 
4.17 0.9324;0.9614 -3.06+85.07i 3.41% 


由 表 1, K 2 可知 ， 四 个 算 例 都 说 明 gOSCR 
对 应 等 效 同 构 系 统 的 最 弱 单 馈 入 系统 与 原 系统 主 
导 特 征 根 轨迹 的 误差 都 较 小 ， 相 对 误差 最 大 为 
0.76%， 是 满足 工程 精度 需要 的 ;而 gSCR 对 应 最 
弱 等 效 单 馈 入 系统 与 原 系统 主导 特征 根 轨迹 的 误 


变 流 器 控制 参数 和 网 络 参 数 见 附录 A 表 A 1。 本 节 
基于 八 馈 入 系统 ， 验 证 等 效 同 构 系 统 和 单机 
COSCR 在 “ 同 构 系统 ”中 的 有 效 性 。 其 中 ， 
VSC1~VSC4 节点 电压 相近 ，VSC5~VSC8 节点 电 
压 相 近 。 因 此 ， 为 节省 表格 空间 ， 后 文 算 例 中 的 节 
点 电压 分 别 为 两 个 代表 节点 电压 值 , 以 反映 节点 电 
压 间 的 差异 。 


一 § 一 


差 较 大 ， 最 大 相对 误差 为 4.96%， 会 对 系统 稳定 性 
分 析 造 成 一 定 误差 。 因 此 ， 在 系统 实际 运行 阶段 ， 
利用 gOSCR 对 应 等 效 同 构 系统 的 最 弱 单 馈 入 系统 
的 主导 特征 轨迹 近似 原 系 统 的 主导 特征 轨迹 是 有 
效 的 ， 用 gSCR 则 无 法 准确 度量 系统 电网 强度 。 
3.2.2 ”临界 广义 运行 短路 比 解析 计算 

为 了 分 析 解 析 计 算 CgOSCR 的 有 效 性 ， 考 虑 
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如 下 算 例 情况 。 

改变 外 部 电网 电压 , 并 增加 线路 电抗 Line910, 
直至 八 馈 入 系统 临界 稳定 , 计算 不 同 电压 下 系统 的 
gOSCR 〈 记 为 “CgOSCR 真实 值 ”); 同时 ， 根 据 
式 计 算 相 应 CgOSCR 〔〈 记 为 “CgOSCR 解析 值 ”)。 
图 4 描绘 了 八 馈 入 系统 运行 在 端口 电压 为 0.98、 
1.05pu 的 临界 稳定 下 , 八 台 电力 电子 设备 端口 电压 
的 等 幅 振 荡 曲 线 。 


1 1.5 2 2.5 3 3.5 4 
t/s 


4 ” 八 馈 入 系统 节点 电压 振荡 曲线 
Fig.4 Terminal voltage oscillation curse of MIPES with 
eight power electronic devices. 
表 3 系统 CgOSCR 真实 值 、 解 析 值 
Table 3 True value and calculated value of CgOSCR 


Ulpu CgOSCR 真实 值 ”CgOSCR 解析 值 ”相对 误差 
0.92;0.99 2.01 2.02 0.50% 
0.98;1.05 2.07 2.03 1.46% 
1.02;1.09 2.06 2.04 0.97% 


# 3 给 出 不 同 临界 电压 下 CgOSCR 真实 值 和 
解析 值 结果 。 可 以 看 出 , 最 大 出 力 及 外 部 电网 电压 
变化 时 ， 系 统 CgOSCR 真实 值 与 解析 值 的 误差 均 
BUN, 最 大 相对 误差 仅 为 1.46%,， 满足 工程 精度 需 
要 ， 验 证 了 解析 计算 CgOSCR 方法 的 有 效 性 。 
3.3 ”哈密 等 比例 缩小 系统 电磁 仿真 分 析 

Cloudpss 是 一 个 基于 云 计算 的 开放 平台 , 用 于 
各 种 功能 , 包括 能 源 互联 网 全 过 程 建 模 和 仿真 、 异 
构 并 行 计算 加 速 等 20P5] 。 

为 利用 电磁 仿真 模型 验证 仿真 计算 系统 
CgOSCR 方法 和 等 效 单 机 系统 的 COSCR 在 “ 同 构 
系统 ”中 的 有 效 性 ， 在 Cloudpss 平台 中 搭建 哈密 
等 比例 缩小 同 构 系统 (网 络 拓 扑 与 实际 哈密 系统 相 
同 , 变 流 器 控制 结构 、 参数 、 容 量 和 网 络 参 数 不 同 ) 
和 相应 等 效 同 构 系统 的 单 馈 入 系统 ， 为 方便 分 析 ， 
上 述 系 统 折算 至 220kV 电压 等 级 。 系 统 拓 扑 如 附 


KA 图 A3(c) 所 示 , 每 片区 中 地 名 下 标 数 字 代表 该 
地 区 中 等 效 变 流 器 编号 ; 其 中 ， 节 点 1-54 为 变 流 
器 设备 馈 入 节点 ,节点 55-93 为 无 源 中 间 节 点 ， 节 
点 94 等 效 为 无 穷 大 电源 节点 ; 变 流 器 参数 和 网 络 
参数 见 附录 A 表 A2。 

调节 线路 电感 Cireoow 和 外 电网 电压 使 哈密 等 
比例 缩小 同 构 系 统 在 不 同 端口 电压 达到 临界 稳定 ， 
计算 不 同 电压 下 系统 的 gOSCR 〈 记 为 “CgOSCR 
AX”) 和 gSCR( 记 为 “CgSCR AXW”); E 
时 , 对 应 调节 等 效 同 构 系统 , 使 其 在 对 应 等 效 端口 
电压 下 达到 临界 稳定 ， 计 算 其 gOSCR Cid W 
“CgOSCR 仿真 值 ”)。 

下 面 进一步 以 哈密 等 比例 缩小 同 构 系统 如 何 
在 额定 工 况 下 达到 临界 稳定 状态 为 例 , 描述 算 例 调 
节 过 程 。 首 先 系 统 稳定 运行 ,在 0.25s 时 ， 对 外 电 
网 施加 -0.05pu 的 电压 跌落 ， 系 统 在 等 效 端口 电压 
为 1 时 达到 了 临界 稳定 状态 。 图 5(a) 描 绘 了 该 过 
程 中 系统 中 第 17 台电 力 电子 设备 端口 三 相 电 压 的 
振荡 波形 ; 图 5 人 b) 描 绘 了 对 端口 电压 进行 快速 全 
里 叶 变 换 的 分 析 , 可 以 看 出 , 系统 振荡 频率 大 致 在 
50+26Hz AA; 图 5(c)、(d) 分 别 描 绘 了 系统 中 所 
有 电力 电子 设备 输出 有 功 功率 、 无 功 功率 的 振荡 波 
形 。 
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5 ”哈密 等 比例 缩小 同 构 系 统 电压 、 有 功 功率 、 无 功 功 
率 振荡 波形 、FFT 分 析 
Fig.5 Voltage, active power, reactive power oscillation 
waveform and FFT analysis of Hami scales down 
homogeneous systems. 
表 4 CgOSCR 真实 值 、 仿 真 值 
Table 4 True value and simulation calculated value of 


CgOSCR 


Uegq/pu 仿真 值 真实 值 相对 误差 


= 


0.9022 1.6320 1.6800 2.86% 
0.9525 1.7298 1.7625 1.86% 
1.0051 1.8312 1.8663 1.88% 
1.0559 1.9239 1.9560 1.64% 
1.0977 2.0150 2.0443 1.43% 

表 5 CgSCR 真实 值 
Table5 True value of CgSCR 
Ueq/pu 0.90 0.96 1.01 1.06 1.10 
真实 值 | 2.0924 1.9696 1.8727 1.7786 1.7202 
误差 17.79% 


K 4 给 出 不 同 临界 电压 下 CgOSCR 真实 值 和 
仿真 值 结果 。 可 以 看 出 ， 外 部 电网 电压 变化 时 ， 系 
统 CgOSCR 真实 值 与 仿真 值 的 误差 均 较 小 ， 最 大 
相对 误差 仅 为 2.86%， 满 足 工 程 精度 需要 ， 验 证 了 
仿真 计算 CgOSCR 方法 的 有 效 性 ， 以 及 利用 单 馈 
入 系统 的 临界 运行 短路 比 能 够 近似 刻画 “ 同 构 系统 ” 
的 运行 稳定 边界 。 
表 5 给 出 不 同 临界 电压 下 CgSCR 真实 值 及 其 
最 大 相对 误差 (不 同 临界 电压 下 CgSCR 真实 值 之 
间 的 最 大 误差 除 以 等 效 端口 电压 为 1.01 的 CgSCR 
真实 值 )。 可 以 看 出 ， 外 部 电网 电压 变化 时 ， 系 统 
CgSCR 真实 值 随 之 发 生 较 大 变化 ， 最 大 相对 误差 
为 17.79%， 不 满足 工程 精度 需要 ， 说 明 在 系统 实 
际 运行 阶段 ， 利 用 CgSCR 刻画 多 馈 入 系统 的 稳定 
边界 存在 较 大 误差 。 
综 上 所 述 , 时 域 、 频 域 仿真 结果 和 广义 运行 短 
路 比 的 解析 分 析 结 果 是 相 一 致 的 , 验证 了 广义 运行 
短路 比 指标 及 其 临界 值 的 有 效 性 
4 结论 


为 考虑 系统 非 额定 运行 方式 (实际 出 力 偏离 额 
定 功率 以 及 外 部 电网 电压 偏离 额定 值 ) 对 系统 电网 
强度 的 影响 ， 本 文 近似 解 厢 了 “ 同 构 系统 ”中 设备 
动态 特性 与 端口 电气 量 的 相互 作用 关系 , 通过 构造 
等 效 同 构 系 统 以 逼近 原 系统 的 稳定 性 ,提出 了 广义 
运行 短路 比 和 临界 广义 运行 短路 比 指标 及 其 实用 
化 算法 ,两 者 的 距离 能 够 更 准确 地 量化 非 额定 运行 
方式 下 系统 的 电网 强度 和 稳定 裕 度 。 广 义 运行 短路 
比 指标 同样 适用 于 度量 “ 异 构 系统 ”的 电网 强度 和 
稳定 裕 度 ， 分 析 方 法 与 “ 同 构 系统 ”类 似 。 
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Generalized Short Circuit Ratio of New Energy Power System Considering Operating State: Generalized 
Operation Short Circuit Ratio 


LIU Chenxi!, XIN Huanhai!, ZHOU Yuhan!, WANG Guanzhong!, YUAN Hui! 
(1. College of Electrical Engineering, Zhejiang University, Hangzhou 310027, China) 


Abstract: The new energy multiple-infeed power electronic system has the issues of sub-/super-synchronous 
oscillation caused by the mismatch between the phase-locked loop bandwidth and the grid strength. It is 
necessary to accurately quantify the grid strength to ensure the stable operation of the system. To this end, this 
paper considers the impact of active power output or terminal voltage deviation from the rated value on the grid 
strength of the system in the operating scenario. Meanwhile, generalized operation short circuit ratio index is 
proposed under non-rated operating conditions. This paper firstly analyzes the change law of the terminal 
characteristics of the devices when the operating mode of the new energy system changes, and approximately 
decouples the interaction between dynamic characteristics of devices and terminal electrical quantities. 
Secondly, the decoupling property of the network-side admittance matrix is studied when considering the 
electrical quantities changing. Considering the operation mode, the generalized operation short circuit ratio 
index and the corresponding practical calculation method are proposed, which expands the application range of 
the generalized short circuit ratio theory. Finally, simulation results verify the effectiveness of the proposed 
index 


Key words: sub-/super-synchronous oscillation; generalized operation short circuit ratio; non-rated operating 
conditions; decoupling property 
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图 A 1 并 网 变 流 器 框图 


Fig A1 Diagram of three-phase grid-connected converter 


图 A2 系统 稳定 裕 度 度量 流程 图 


4J 
Fig.A2 Flow chart of system stability margin measurement 
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A3 单 馈 入 、 八 馈 入 和 哈密 等 比例 缩小 系统 仿真 模型 
Fig.A3 The test system with the VSC , eight VSCs and Hami scales down homogeneous systems. 


RAL 时 域 仿真 模型 参数 


Table A1 Parameters of time domain simulation 


变 流 器 控制 参数 
SukVA 1500 Ud/V 690 Lepu 0.05 Cepu 0.05 
Kip 0.2 Ki 10 了 Piewpu 1 Qree/pu 0 
Tr 0.0001 Kopiip! 12 Koni 7200 
网 络 参数 
Line 12/pu 0.1 Line 23/pu 0.05 Line 34/pu 0.035 Line 49/pu 0.02 
Line 56/pu 0.2 Line 67/pu 0.15 Line 78/pu 0.018 Line 89/pu 0.135 


表 A2 Cloudpss 电磁 仿真 模型 参数 


Table A2 Parameters of electromagnetic simulation model on Cloudpss 


变 流 器 控制 参数 


SukVA 1000 UV 620 Li/pu 0.1225 Cr/pu 0.0604 
Kacp 5 Kaci 20 Kip 0.5 Ki 50 
Pres/pu 1 Qret/pu 0 Koiip2 25 Konia 18000 

Cac/pu 0.1089 Uaco/V 1100 

网 络 参数 

Linelss/pu 0.067 Line2,ss/pu 0.067 Lines,ss/pu 0.067 Line4,ss/pu 0.067 
Lines.s6/pu 0.067 Lineo,s6/pu 0.067 Line7.s7/pu 0.067 Lines,ss/pu 0.067 
Lines,ss/pu 0.067 Line10,59/pu 0.067 Line11,60/pu 0.067 Line12,60/pu 0.067 
Line13,61/pu 0.067 Linel4.62/pu 0.067 Line1s,62/pu 0.067 Line 16,63/pu 0.067 
Line17,64/pu 0.067 Line1s,65/pu 0.067 Line19,65/pu 0.067 Line20,66/pu 0.067 
Line?21,67/pu 0.067 Line22,68/pu 0.067 Linez3.68/pu 0.067 Line24,68/pu 0.067 
Line25,68/pu 0.067 Line26,69/pu 0.067 Line27,69/pu 0.067 Line2s,70/pu 0.067 
Line29,70/pu 0.067 Lineso,71/pu 0.067 Linesi,71/pu 0.067 Line32,72/pu 0.067 
Line33,73/pu 0.067 Lines4,73/pu 0.067 Liness,74/pu 0.067 Line36,7s/pu 0.067 
Lines7,76/pu 0.067 Linesg,77/pu 0.067 Line39,7s/pu 0.067 Line4o,7s/pu 0.067 
Linea1.78/pu 0.067 Line42,78/pu 0.067 Line43,79/pu 0.067 Linea4.so/pu 0.067 
Lineas.so/pu 0.067 Line46,80/pu 0.067 Line47,80/pu 0.067 Line4s,81/pu 0.067 
Line49s2/pu 0.067 Lineso,82/pu 0.067 Lines1,83/pu 0.067 Lines2,84/pu 0.067 
Lines3,84/pu 0.067 Lines4,ss/pu 0.067 Liness,86/pu 0.00018282 Lines6,86/pu 0.00076013 
Lines7,86/pu 0.0017608 Liness,86/pu 0.0024632 Lineso,86/pu 0.0021264 Line6o,s6/pu 0.0014818 
Linesı,86/pu 0.0012701 Line62,86/pu 0.00066391 Linees.87/pu 0.000789 Line6487/pu 0.0023381 
Liness.s7/pu 0.0017416 Linees.s7/pu 0.0014048 Linee7.g7/pu 0.0010777 Linees.s7/pu 0.00025979 
Lineeo.s7/pu 0.00011546 Line70,87/pu 0.00037526 Line7lsspu 0.0018948 Line72,8s/pu 0.0016357 
Line73,38/pu 0.0011398 Line74,s8/pu 0.000977 Line7s.ss/pu 0.00060692 Line76.ss/pu 0.000977 
Line77.8s/pu 0.00060692 Line7s,90/pu 0.00014063 Line79.90/pu 0.0016357 Lineso,89/pu 0.00014063 
Linesı,89/pu 0.0016357 Lines2.91/pu 0.0011398 Lineg3,91/pu 0.000977 Linesa,o1/pu 0.0011398 
Linegs,91/pu 0.000977 Linege,s7/pu 0.005181 Linego,90/pu 0.005181 Lines7,92/pu 0.007789 | 
Linegs,92/pu 0.0032454 Lineso,92/pu 0.0058418 Lineo1,92pu 0.00064909 Lines2,93/pu 0.0025173 
Lineo3,94/pu 0.003 


